Evolution from BCS to BEC superfluidity in the presence of spin-orbit coupling 
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We discuss the evolution from BCS to BEC superfluids in the presence of spin-orbit coupling, and 
show that this evolution is just a crossover in the balanced case. The dependence of several thermo- 
dynamic properties, such as the chemical potential, order parameter, pressure, entropy, isothermal 
compressibility and spin susceptibility tensor on the spin-orbit coupling and interaction parameter 
at low temperatures are analyzed. We studied both the case of equal Rashba and Dresselhaus (ERD) 
and the Rashba-only (RO) spin-orbit coupling. Comparisons between the two cases reveal several 
striking differences in the corresponding thermodynamic quantities. Finally we propose measuring 
the spin susceptibility as a means to detect the spin-orbit coupling effect. 

PACS numbers: 03.75.Ss, 67.85.Lm, 67.85.-d 



Superfluidity is a ubiquitous phenomenon that is en- 
countered in nearly every area of physics including con- 
densed matter, nuclear, astro, and atomic and molec- 
ular physics. Superflow results from strong correla- 
tions between particles, which for any given interacting 
Fermi system could not be controlled externally until re- 
cently with the advent of ultra-cold atoms. In standard 
condensed matter there is a continuous search for new 
charged superfluids (superconductors) since the type and 
strength of interactions can not be tuned even within the 
same class of materials. In the case of nuclear matter 
the issue of tunability of interactions is even worse, be- 
ing hopeless for neutron stars. However, the situation is 
much more favorable for ultra-cold Fermi atoms, where 
the ability to control interactions between particles, via 
Feshbach resonances, has been demonstrated in exper- 
imental studies of the so-called crossover from BCS to 
BEC superfluidity 

Further control of interactions is now possible through 
newly developed experimental techniques that allow the 
production of fictitious magnetic fields which couple to 
neutral bosonic atoms [3, 0] ■ These fictitious magnetic 
fields are generated through an all optical process, but 
produce real effects like the creation of vortices in the 
supcrfluid state of bosons. Furthermore, artificial spin- 
orbit coupling has also been produced in neutral bosonic 
systems [3| where the strength of the coupling can be con- 
trolled optically. In principle the same techniques can be 
applied to ultracold fermions 0, H[ , which, when coupled 
with the control over the interaction using Feshbach res- 
onances, allows for the exploration of superfluidity not 
only as a function of interactions, but also as a function 
of fictitious magnetic fields [H[ , or as a function of spin- 
orbit coupling discussed here. An introduction to the 
effects of controllable fictitious magnetic and spin-orbit 
fields can now be found in the literature @. 

It is in anticipation of experiments involving spin-orbit 
coupling in fcrmionic atoms such as 6 Li, 40 K and isotopes 
of Ytterbium, that we discuss here the evolution from 
BCS to BEC superfluidity in the presence of controllable 



spin-orbit couplings for balanced fermions in three di- 
mensions. We investigate spin-orbit effects with Dressel- 
haus and/or Rashba 0] terms, and analyze several 
thermodynamic quantities including the order param- 
eter, chemical potential, thermodynamic potential, en- 
tropy, pressure, isothermal compressibility, and spin sus- 
ceptibility tensor as a function of spin-orbit coupling and 
interaction parameter at low temperatures. We conclude 
that the BCS-to-BEC evolution for balanced fermions in- 
cluding spin-orbit effects is just a crossover. 

Hamiltonian: To address the problem of the evolu- 
tion from BCS to BEC superfluidity in the presence of 
spin-orbit fields for balanced or imbalanccd Fermi-Fermi 
mixtures, we start with the generic Hamiltonian density 



H(r) = H (r)+Hi(r). 
The single-particle Hamiltonian density is 



(1) 
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where K a — — V 2 /(2m Q ) — \i a is the kinetic energy in ref- 
erence to the chemical potential [i a , and hi (r) is the spin- 
orbit field along the i-direction (a =f, = x, y, z). The 
interaction term is Hi(r) = — g , 0|(r)'0j(r)^(r)i/)-f (r), 
where g is a contact interaction. In this paper we set 
h = ks = 1- 

Effective Action: The partition function at tempera- 
ture T is Z = J T>[ip,ipi]faq? (-S[ip, with action 



= I drdr 







(3) 



Using the standard Hubbard-Stratanovich transfor- 
mation that introduces the pairing field A(r,r) = 
g(ipl(v, T)xj}f(r, t)) we can write the intermediate action 
Stoh^A.At] = S ro [V>,V t ] + SihM t ,A,At], where 
the no-interaction action is 



S'no[V' J ' ! / ,t ] 



drdr 
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and the action due to the auxiliary field is 
|A(r,r)| 2 



Si 



drdr 
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Using the four-dimensional vector ' 
{^jj V'ti V'ili the intermediate action becomes 



drdr 



!A(r,r)| 
.9 



The 4x4 matrix M is 

( d T + A' t —h± 

-hi d T + k x 

At 
\ -At 



M = 
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where h± = h x — ih y corresponds to the transverse com- 
ponent of the spin-orbit field, h z to the parallel com- 
ponent with respect to the quantization axis z, = 
— hz, and K 



which describes the Hamiltonian of the elementary ex- 
citations in the four-dimensional vector basis 'ft = 
|^|(k),^l(k),-0f(— k),-04,(— k)} . The spin-orbit field is 

hj_(k) = hzj(k) + h£j(k), where the first term is of the 
Rashba-type h#(k) = Vr (— k y i + k x y) , and the sec- 
ond is of the Dresselhaus-type ho(k) = (k y ic + k x y) . 
We assume, without loss of generality, that vr > and 
Vr> > 0. The magnitude of the transverse field is then 

hj_(k) = J (v D - v R ) 2 k 2 , + (v D + v R ) 2 k%. In the lim- 
iting cases of pure Rashba (R) with Vd = and for 
equal Rashba-Drcssclhaus (ERD) couplings with vr = 

vd = v/2, the transverse fields are hj_(k) = VrJ~K%. + k 2 . 

(vr > 0) and h±(k) = v\k x \ (v > 0), respectively. 

Order parameter and number equations: The saddle 
point thermodynamic potential £Iq = TSq is obtained by 
integrating out the fcrmions leading to 

T 



^> and ^t leads to the effective action 
|A(r,r)| 2 



Integration over the fields fl Q = V— ^ In {1 + exp [-Ej(k)/T]} + ^ K + , 
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J [ g 2V T 

where K+ = (A t + 

Saddle Point Approximation: To proceed we use the 
saddle point approximation A(r, r) = A + 7](t,t), and 
separate the matrix M into two parts. The first one is 
the saddle point matrix Mo, where the transformation 
A(r, t) — > Aq takes M — > Mo- The second one is the 
fluctuation matrix Mp = M — Mo, which depends only 
on r)(r, r) and its Hcrmitian conjugate. 

Using the saddle point approach we write the effective 
action as Seff = So + Sf, where 

|2 



with K + = K^(— k) + K_i(—k) /2. The order parame- 
ter is determined via the minimization of fig with respect 
to |A | 2 leading to 



V 1^ 
-=--^n F [E 3 .(k} 



dE 3 (k) 
WoJ'' 
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Sq = drdr 
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where n F [E 3 (k)\ = l/(exp [E 3 (k)/T] + 1) is the Fermi 
function for energy Ej(k). We replace the contact in- 
teraction g by the scattering length a s through the rela- 
tion l/g = -m + /(47ta s ) + (1/V)£k [l/(2ek,+)] , where 
m + = 2m±m^/(m\. + m^) is twice of the reduced mass, 
£k.Q = k 2 /(2m a ) are the kinetic energies, and £k,+ = 
[ £ k,t + e k.jJ /2. The number of particles at the saddle 
point is obtained by N a — —d£lo/d/j a , leading to 



is the saddle point action and 
>?(r,r)| 2 T 



St 



drdr 



In dot ( 1 



M„ 1 M f 



9 2U 

is the fluctuation action in all orders in the fluctuation 
field. The effects of fluctuations at both zero temperature 
and near the critical temperature will be discussed later. 

A transformation to the momentum-frequency coordi- 
nates (k,iu) n ), where u„ = (2n + 1)tvT, leads to 



(8) 
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The self-consistent relations shown in Eqs. (J7J and © 
are general for arbitrary mass and population imbalances. 
However, next, we particularize our discussion to the case 
of a balanced system with equal masses. 

Balanced Populations: In the case of mass and pop- 
ulation balanced systems, the four eigenvalues of the 

matrix Hq are £a(k) = 



Mk)\ 



E 2 {k) 



where Ej (k) are the eigenvalues of the matrix 



Hn 



( ^t(k) 
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e 2 (kr + |A | 2 , E 3 (k) = -E x (k), and E 4 (k) = 
—E2(k). Here, the auxiliary energies are £i(k) = ^(k) + 
hj_(k), and £2(k) = ^(k) — h±(k). The corresponding 
order parameter equations at the saddle point level is 



(0) 



7 ~ 2 ^ 

y k 



Xi(k) X 2 (k) 



2E x {k) 2E 2 {k) 



(9) 
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where X m (k) = tanh [E m (k)/2T] (m = 1,2). Since the 
mixture of equal mass fermions is balanced, the chemical 
potentials are the same m = = /i, and the associated 
number equation is N = —dQ/dfi that reduces to 

*=E[i-i§^>-t§-4 do) 

In Fig. [TJ we show the zero temperature behavior of 
Ao| and ^ as a function of l/(fcpa s ) for various values 
of spin-orbit coupling in the equal-Rashba-Dresselhaus 
(ERD) and for Rashba-only (RO) cases. In the ERD case 
the order parameter |A | is independent of v, and the 
chemical potential n(v) is simply (i(v) = (i(0) — mv 2 /2, 
since the transverse field hj_(k) = v\k x \ can be elimi- 
nated by momentum shifts along the ^-direction, effec- 
tively gauging away spin-orbit effects in the charge or 
momentum sector. This symmetry also implies that the 
critical temperature T c as a function of l/(kpa s ) for finite 
v is the same as that for v = 0. However, in the RO case, 
shifts in momentum can not gauge away the spin-orbit 
coupling, and |Ao| increases with increasing vr, while /i 
decreases as vr increases, exhibiting the same tendency 
as in the ERD case. In the BCS regime, the increase of 
Ao| with vr also leads to an increase of T c with increas- 
ing v R . 




l/k F a s l/k F a s 



FIG. 1: Order parameter |Ao| and chemical potential /i (in 
units of the Fermi energy ef) as a function of interaction pa- 
rameter l/(kF<is) for different spin-orbit couplings vr/vf — 
(solid), Vr/vf = 0.8 (dashed), vr/vf = 1.0 (dotted), and 
Vr/vf = 1.2 (dot-dashed) at T — in the RO case. Here 
Vf = kF/m is the Fermi velocity. 

Momentum distribution and excitation spectrum: The 
momentum distribution n(k) is obtained from Eq. (|10[) 
using the definition N = X)k n (^)- At fixed momentum 
component k z = and fixed interaction strength, the 
momentum distribution n(k) shifts continuously with in- 
creasing spin-orbit coupling in the BCS \\/(kpa s ) <C — 1] 
or unitarity regimes [l/{kpa s ) — > 0], For zero spin-orbit 
coupling, n(k) is that of a supcrfluid degenerate Fermi 
system with identical single-particle bands £(k) and has 
a nearly flat momentum distribution until the Fermi mo- 
mentum is reached. However, as the spin-orbit cou- 
pling is turned on, non-identical single-particle bands 
^(k) = £(k) - h ± (k) and &(k) = £(k) + h ± (k) in 
the helicity basis |ff), |JJ.) emerge and produce a double 



structure with a reasonably flat momentum distribution 
centered around finite momenta in the (k x ,k y ) plane. In 
the BEC regime [l/(kpa s ) 3> 1] the momentum distribu- 
tions for weak and strong spin-orbit coupling broadens 
substantially due to the loss of degeneracy in the Fermi 
system when the chemical potential goes below the min- 
ima of the helicity bands and becomes large and nega- 
tive. Even though there is a substantial change in the 
momentum distribution as a function of the spin-orbit 
coupling, we notice that the excitation energies E\(k) 
and E2(k) is always gapped for all values of the inter- 
action parameter l/(kpa s ) or the spin-orbit field h±(k), 
immediately suggesting that thermodynamic properties, 
which depend on the excitation energies, evolve smoothly 
from the BCS to the BEC regime in the balanced case 
for fixed values of spin-orbit coupling. The omnipres- 
ence of a gap in the excitation spectrum shows that the 
evolution from BCS to BEC superfluidity at finite spin- 
orbit coupling for balanced systems is just a crossover. 
The situation is different for imbalanced systems, where 
gapless regions emerge in the excitation spectrum and 
topological phase transitions occur, so long as the sys- 
tem is stable 0, [l(j ■ A thermodynamic signature of this 
crossover for balanced systems is seen in the isothermal 
compressibility discussed next. 

Isothermal compressibility: An important thermody- 
namic property, which can now be measured experimen- 
tally using the fluctuation-dissipation theorem, is the 
isothermal compressibility 

1 fdP\ V fdN\ 

As shown in Fig. [2^,, for the RO case, the isother- 
mal compressibility kt at fixed interaction parameter 
l/(kpa s ) increases with increasing spin-orbit coupling 
vr, as the Fermi system becomes less degenerate reducing 
the Pauli pressure, and thus more compressible. How- 
ever, in the ERD case, the isothermal compressibility for 
fixed interaction parameter does not change with increas- 
ing spin-orbit coupling v. In this high symmetry situa- 
tion the momentum shift in the energy spectrum and the 
accompanied shift in the chemical potential do not affect 
the degeneracy of the Fermi system or the Pauli pres- 
sure, leading to an isothermal compressibility which is 
independent of the spin-orbit coupling v. 

Equation of State and Entropy: Since the thermody- 
namic potential 51 = — PV, the saddle point pressure 
is Po(T,/j, a ) = — Clo/V, which can be shown to be al- 
ways positive for arbitrary spin-orbit coupling. The gen- 
eral trend of the pressure for fixed interaction param- 
eter (from the BCS to the unitarity regimes) is to de- 
crease with increasing spin-orbit coupling for both ERD 
and RO cases. The situation in the BEC regime re- 
quires the inclusion of quantum fluctuations to recover 
the corresponding Lcc-Yang corrections in the presence 
of spin-orbit effects. The entropy is then calculated from 
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FIG. 2: a) Compressibility kt (in units of l/(neF)) as a 
function of interaction parameter l/(fc_Fa 3 ) at T = and b) 
entropy per particle So /N as a function of temperature T (in 
units of ep) at unitarity in the RO case, for the same group 
of spin-orbit coupling values vr/vf as in Fig. [T] 



S = — (8Q/dT) v . In Fig.[2p, we show the saddle point 
entropy So for the RO case at unitarity. For fixed T, So 
decreases with increasing spin-orbit coupling due to the 
stabilization of superfluidity by the spin-orbit field. 

Spin Susceptibility Tensor: A rotation of the matrix 
Ho into the helicity basis |fr), j-JJ.) introduces order pa- 
rameters Ao,fff and Ao,^^,, which are controlled by the 
spin-orbit coupling. The emergence of the triplet com- 
ponent affects dramatically the spin susceptibility of the 
system. Using standard linear response theory [11|, the 
uniform spin susceptibility tensor per unit volume is 



Xij 



2 

Ms \ - 



V 



k 



[oy(k)-&y(k)] 



(12) 



where the spin-spin correlations in the single-particle 
channel are ajj(k) = JY W Tr [<TjG(k, iu>)ajG(k, ioj)] 
and in the pair (anomalous) channel are bij(k) = 
E^Tr [cr 4 F(k,ia;)crjFt(k,^)] . The matrices G and F 
are the block matrices appearing in the inverse of M de- 
fined in Eq. (g]), 



i.e., x(T — > 0) = 0. The existence of non-zero spin re- 
sponse even at T = is a direct measure of the induced 
triplet component of the order parameter due to the pres- 
ence of spin-orbit coupling, since that a pure singlet su- 
pcrfluid at T — must have zero spin susceptibility since 
all fermions are paired into a zero-spin state. 




FIG. 3: a) Spin susceptibility \ zz (in units of n B n/tF) as a 
function of l/(fcFa s ) at T = in the ERD case for v/vf = 
(solid), v/vf = 0.8 (dashed), v/vf = 1 (dotted), and v/vf = 
1.2 (dot-dashed), b) Spin susceptibility \ zz as a function of 
v/vf at T = at unitarity in the ERD case. 



Conclusions: We have studied the effects of spin-orbit 
coupling in the evolution from BCS to BEC superfluidity 
at low temperatures, and concluded that this evolution 
is just a crossover. We discussed effects of spin-orbit cou- 
pling on thermodynamic properties including the order 
parameter, chemical potential, pressure, entropy, isother- 
mal compressibility and spin susceptibility tensor to sup- 
port the crossover picture. We also proposed way to 
experimentally detect the spin-orbit coupling effect by 
measuring the spin susceptibility. 
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In Fig. [3K, we show plots of Xzz for the ERD case at 
T = as a function of l/(kpa s ) for various values of spin- 
orbit coupling, and the behavior of Xzz for the RO case 
is qualitatively similar. In Fig. |3p, we show Xzz versus v 
in the unitary limit l/(kpa s ) = 0. The maximum in \ zz 
corresponds to the maximum in the triplet component of 
Ao. For small and large v the triplet component is small. 

In the ERD case \zz = Xxx ^ Xyy> an d in the zero 
temperature limit Xyy{T — > 0) = 0, while \zz = Xxx 
remains finite for non-zero spin-orbit coupling. In the 
RO case Xzz ^ Xxx = Xyy> an d m the T — > limit 
Xxx(T 0) = Xyy (T -)• 0) = X zz(T -> 0)/2. Lastly, 
for hj_(k) = (no spin-orbit coupling) the spin sus- 



ccptibilty tensor becomes Xij 
X = [fi 2 B /(2VT)]J2 kS ech 2 



X&ij-, where the scalar 



is the 



|A | 2 /(2T) 



Yoshida function, which vanishes at zero temperature, 
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